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I. INTRODUCTION 



Bulk inflaton shadows of vacuum gravity 

Tsutomu KobayashQ and Takahiro Tanak£0 
Department of Physics, Kyoto University, Kyoto 606-8502, Japan 

We introduce a (5 + m)-dimensional vacuum description of five-dimensional bulk inflaton models 
with exponential potentials that makes analysis of cosmological perturbations simple and transpar- 
ent. We show that various solutions, including the power-law inflation model recently discovered by 
Koyama and Takahashi, are generated from known (5 + m)-dimensional vacuum solutions of pure 
gravity. We derive master equations for all types of perturbations, and each of them becomes a sec- 
ond order differential equation for one master variable supplemented by simple boundary conditions 
on the brane. One exception is the case for massive modes of scalar perturbations. In this case, 
there are two independent degrees of freedom, and in general it is difficult to disentangle them into 
two separate sectors. 
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\£) • Recent progress in particle physics suggests that the universe might be a four-dimensional subspace, called a "brane" , 
embedded in a higher dimensional "bulk" spacetime. In this braneworld picture, ordinary matter fields are supposed 
to be confined to the brane, while gravity can propagate in the bulk. Various kinds of braneworld models have been 
proposed, and the cosmological consequences of these models have been studied (for a review see, e.g., Ref. [l]). The 
idea of the braneworld brings new possibilities, in particular to scenarios of the early universe. 

A simple model proposed by Randall and Sundrum 0, is such that the unperturbed bulk is a five-dimensional 
anti-de Sitter spacetime (warped bulk) bounded by one brane or two. A homogeneous and isotropic cosmological 
solution based on this model has been explored 0, la, IE S ■ A slow- roll inflation driven by a scalar field confined to 
' the brane was considered in An empty bulk, however, seems less likely from the point of view of unified theories, 
which often require various fields in addition to gravity. Considering a bulk scalar field, Himemoto et al 0, 0, IT^ 
have shown that, interestingly, a bulk scalar field can mimic the standard slow-roll inflation on the brane under a 
! certain condition (see also Ref. 0]). 

1 Also in the context of heterotic M theory cosmological solutions has been studied 0, 0, 0, . In the model 
Q_i! discussed in Refs. 0, OH, the scalar field has an exponential potential in the bulk and the tensions of the two 
branes are also exponential functions of the scalar field. In this model the power-law expansion (but not inflation) is 
realized on the brane. A single-brane model with such exponential-type potentials is also interesting, and it has been 
investigated for a static brane case 0, 0, |2(J and a dynamical (cosmological) ca se |2ll | 22l I2H. l24j . Very recently, 
an inflationary solution was found in a similar setup by Koyama and Takahashi |25l l2fj||. extending the results of 
Refs. |l9ll2Ct l2lLl22| . A striking feature of their model is that cosmological perturbations can be solved analytically. 

In spite of the tremendous efforts by many authors (23, |28|, |2|| |3(J H3, H3, H3, HJ] , it is still an unsolved problem 
to calculate the evolution of perturbations in the braneworld models with infinite extra dimensions. This lack of 



knowledge constrains the predictability of this interesting class of models. There is an approximate way to estimate 
the density fluctuations evaluated on the brane, in which perturbations in the bulk are neglected [jj. However, once 
we take into account perturbations in the bulk, generally we cannot avoid solving partial differential equations in the 
bulk with discouragingly complicated boundary conditions. 

Only a few cases are known where perturbation equations can be analytically solved [32l l33l l34| . One of them 
is the special class of bulk inflaton models mentioned above (2j| [2(| . In this paper we clarify the reason why the 
perturbation equations are soluble in this special case. Based on this notion, we present a new systematic method to 
find a wider class of background cosmological solutions and to analyze perturbations from them. 

This paper is organized as follows. In the next section, we explain our basic ideas of constructing background 
solutions and of analyzing cosmological perturbations. In Sec. Ill we consider a model with a single scalar field in the 
bulk, which is the main interest of this work, and derive an effective theory on the brane. Then, in Sec. IV we present 
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some examples of exact solutions for the background cosmology obtained by making use of the ideas explained in 
Sec. II. Section V deals with cosmological perturbations. Section VI is devoted to discussion. 



II. BASIC IDEAS 
A. Solutions in the Randall- Sundrum vacuum braneworld 

We begin with a model whose action is given by 

S = S g + S b , (2.1) 

where 

S g = -4— [ d D+1 xV^G (R[G] - 2A D+1 ) , (2.2) 

^ K D+l J 

is the action of (D + l)-dimensional Einstein gravity with a negative cosmological constant A£> +1 = —D(D — l)/2£ 2 , 

S b = - [ d D X^J~g~ a, (2.3) 



is the action of a vacuum brane with a tension er, and g is the determinant of the induced metric on the brane. 
We assume Z2 symmetry across the brane, so that the tension of the brane is determined by the junction condition 

as 

2 ^ D I 1 ^* 

where Hq is related to the Z?-dimensional cosmological constant induced on the brane by 

k b = \{D-l){D-2)Hl (2.5) 

and it represents the deviation of a from the fine-tuned Randall- Sundrum value, 2{D — 1) / ln 2 D+1 . 

One of the key ideas in the present paper is to make use of the following well-known fact. If a metric ds 2 D ^ is a 
solution of the D-dimensional vacuum Einstein equations with a cosmological constant At,, 

dsf D+1) = e 2 ^(dz 2 + dsf D) ), (2.6) 

is a solution of the (D + l)-dimensional model defined by Eq. I|2.1|l . and the warp factor is given by 

(2.7) 



smh(Hoz) 

(For a Ricci-fiat brane, we have A& = 0. In this case the warp factor reduces to e w(z ) = l/z.) Namely, we can 
construct a (D + l)-dimensional solution in the Randall-Sundrum braneworld from a vacuum solution of the D- 
dimensional Einstein equations. A well-known example is the five-dimensional black string solution obtained from the 
four-dimensional Schwarzschild solution 1.351. 



B. Bulk inflaton models from dimensional reduction 



We explain how to obtain an (n + 2)-dimensional braneworld model with bulk scalar fields from (n + 2 + ^ ji)[= 
D + l]-dimensional spacetime by dimensional reduction. We use n to represent the number of uncompactified spatial 
dimensions on the brane, which is three in realistic models. Let us consider (n + 2 + ^2 ji)-dimensional spacetime 
whose metric is given by 

ds2 (D+i) = G AB dX A dX B = G ab {x)dx a dx b + e 2Mx) da 2 , (2.8) 
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where daf is the line element of a ji-dimensional constant curvature space with the volume Vj. Here indices a and b 
run from to n + 1, and ^ is assumed to depend only on the (n + 2)-dimensional coordinates x a . 
Then dimensional reduction to n + 2 dimensions yields 



c(n+2) _ ! f m+2 

9 ~ 9k 2 

ZK n+2 J 



xV—G e 



Q 



R[5] -J2^ gabd «& db & + 5 ab d a Qd h Q 



where 



-2A D+1 + ^Kijitii-Ve 



2<f. : 



(2.9) 



(2.10) 



K n+2 := ^ri+i/n^i an d Ki represents the signature of the curvature of the metric daf: —1 (open), (flat) or 1 
(closed). Making a conformal transformation to the "Einstein frame", 



(2.11) 



we have 



o(n+2) _ 

9 ' 2nl 



n+2 J 



R[9] - £>0 o6 0„&&& - -G ab d a Qd b Q 



(2.12) 



Notice that the kinetic term of each scalar field has an appropriate signature since ji > 0. A parallel calculation gives 



-,(71+2) 



x^f — q ae 



-Q/n 



(2.13) 



where a = a\\Vi, so that = K i)+i 'i an d 1 is the determinant of the induced metric on the brane, := 

&ab\ z = z b — e2 ®^ n &ab\ z=z b- The hat upon indices represents restriction to the subspace parallel to the brane. Hence a 
and b run from to n. Here z b represents the location of the brane. In this manner, we can derive (n + 2)-dimensional 
braneworld models with bulk scalar fields which have exponential-type potentials both in the bulk and on the brane. 



III. SINGLE SCALAR FIELD IN THE BULK 



From now on, for simplicity, we focus on models with a single bulk scalar field. The m-dimensional space represented 
by e 2 ^da 2 is compactificd either on a torus (K = 0), on a sphere (K = 1), or on a compact hyperboloid (K = — 1). 
Using a canonically normalized field tp := n~+ 2 \J m(m + njjn <fi, the (n + 2)-dimensional reduced action is written as 

S {n+2) = J d n+2 x\[^ R[Q) - l -G ab d a ^dw - U(^) j - J d n+1 x^—q U(<p), (3.1) 

where the potentials are 

_ (m + n)(m + n + I) 2V z bKn+2V Kmjm - 1) ^ n+ay/n6 , . 

V W> ~ K 2 02 e 9k 2 e ' 

/jK n+2 t /jK n+2 

U(if) = ae-^ 6re "+ 2V , (3.3) 

with 
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FIG. 1: Schematic of the higher dimensional vacuum description of an (n+2)-dimensional bulk inflaton model and the derivation 
of the (n + 1) -dimensional effective theory on the brane. The top-left picture represents the bulk inflaton model which we are 
interested in. To analyze cosmological background solutions (and "zero mode" perturbations), we use the (n+ l)-dimensional 
description shown in the bottom-left corner. On the other hand, we make use of the (D + l)-dimensional description presented 
in the top-right corner to simplify the perturbation analysis. 



If we assume that the metric ds 2 D+1 ^ is given in the form of (|2.fci|) . the action can further reduced to the (n + 1) 
dimensional effective one on the brane. We write the D-dimensional part of the metric in the form of 



ds\ D) = g M dx A dx B = g2 +1 \x)dx a dx b + e 2a{t) do 2 . (3.5) 
Then, comparing the coefficient in front of da 2 , 

<f>(t,z):=a(t)+u(z), (3.6) 
follows. Also, [n + 2)-dimensional metric Q ab is written as 

g ab dx a dx b = e 2m4,/n g ab dx a dx b = e 2 ( m +"W" (q, ah dx h dx h + e 2ma/n dz 2 ^ , (3.7) 

where we have set e"*- 26 -* = 1. Substituting the above expression into l|3.1f) . we perform the integration over z to obtain 

I f , f ™„ /„ „r , m(m + n) 



i n+l x^[e ma/n R[q] 



? ma/n q &b daad b 



-2k b e- ma ' n + Km(m - l) e - 2a - raa /"J ; 



(3.8) 



where 



with 



K n+i — k d/Vi 



\D+1 



e (m+n)u dz 
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and V is the volume of the m-dimensional compactified space. 

The above reduction to (n+1) dimensions can be done more easily, starting with the (D+l)[= n+2+m]-dimcnsional 
action. First we perform the integration over z and obtain a £>-dimensional effective action, 

S {D) = A" / d D x^—g{R[g] - 2A h ), (3.9) 

D J 



where A& is, as before, the one defined by Eq. I|2.5[) . The obtained effective action is that for D-dimensional pure 
gravity with a cosmological constant A^. Compactifying m-dimensions further, and taking into account 

1ai = e 2ma/n 9 ( 2 +1 \ (3-10) 

the same expression for the (ri + l)-dimensional effective action 1)3.8(1 can be recovered in a parallel way as we did for 
the reduction from (D + 1) dimensions to (n + 2) dimensions. Using <*pbd '■= e mQ /™ and ujbd '■= \l n(m + n) /to, we 
can rewrite the action into a familiar form: 

= / d n+1 xj~q L BD R - ^(^ BD ) 2 - 2K Wb I + Km(m l)^ +2 " )/m j , (3.11) 

which implies that the effective theory on the brane is described by a scalar-tensor theory. The action 1)3.8(1 . or 
equivalently (|3.11|) . describes not only the background unperturbed cosmology but also the zero mode perturbation, 
both of which are independent of the extra-dimensional coordinate, z, apart from the overall factor e 2u> . 

The induced metric on the brane q a b is the metric in the Jordan frame. If we use the metric in the Einstein frame, 

dab = e*r a ^, (3.12) 

the effective action becomes 

S in+1) = J d n+1 x^{^±-R[q] - \t h d & (pd^ - V(tp)} > (3-13) 

where <p :— \J m(m + n)/n a and the potential is 

K 2 n+1 Vm = A b e-^»^ ~ Km ^- l ) e - ^X + -r ] ^. (3.14) 

Obviously, the system defined by the above action is equivalent to Einstein gravity with a scalar field. A discussion 
on this type of potential can be found in a recent paper by Neupane [3jj . 

IV. EXAMPLES OF THE BACKGROUND SPACETIME 

In this section, we give some examples of D-dimensional vacuum solutions, which generate (D + l)-dimensional 
braneworld solutions by making use of the prescription described in Sec. Ill Al Here we discuss models with a single 
scalar field and investigate their cosmological evolution in detail. Generalization to the case of multi scalar fields is 
given in Appendix B. 

A. Kasner type solutions 

We first consider the following Kasner-type solution as an example of the Ricci-flat case, At, = 0, 

g Af3 dx A dx § = e 2a{v) [-dif + -y^dy^dy"] + e^^S^dx'dx 3 , 

where 7^ is the metric of m-dimensional constant curvature space, and i and j run from 1 to n. Under the assumption 
of the above metric form, we solve the vacuum Einstein equations, 



e 2a R^ = n 



(5' 2 - a 1 (3 1 + (3" +ma" = 0, (4.1) 
K(m - 1) + noi '0 + (m - l)a' 2 + a" = 0, (4.2) 
2a Rr> = §J n f3' 2 + (m - l)a'(3' + /?"! = 0, (4.3) 



2 a / > v 
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where prime denotes differentiation with respect to 77. Eliminating a" and a' from the above three equations, we have 



n(m +71—1) 



/3' 4 - Kim,- l) 2 /3' 2 . 



We can easily integrate this equation. For example, when K = 1 (compactified on m-sphere S m ), the solution of this 
equation becomes 



n/3' = 



±(m — l)q 
sin[(m — 1)77] 



where 



, + n-V 

and the integration constant was used to shift the origin of time. Thus we find 



(4.4) 



O0 



tan 



f 



±q/n 



Substituting this result into Eq. (|4.3|1 . we have 

a' = cot[(m — l)r]} =p 

This is integrated as 

e (m-i)a = gin [( m _ ^ 



sin[(m — 1)77] 



cot 



m — 1 



-V 



The solution for K = — 1 is easily obtained by replacing sin, tan and cot, in the above expressions by sinh, tanh and 
coth, respectively. The solution for K = behaves like e 13 cx -q ±q l n and e( m_1 ) Q cx r/ 1=F9 . 

Let us further investigate the cosmology of the above example. Setting n — 3 and K = 1, the induced four- 
dimensional metric becomes 



q.~ h dx a dx h 



„2<> 



where 



e" = {sin[(m - 1)77]}^ {cot[(m - 1)77/2]}' 
= {tan[(m- l)77/2]} 9/3 , 



(4.5) 



(4.6) 



and q = y/3m/(m + 2) |45| . We set q to be positive without any loss of generality since the signature of q is flipped 
by a shift of the origin of time, 7/^7/ + n. Here one remark is in order. In the original (D + l)-dimensional model 
m represents the number of the compactified dimensions and therefore is supposed to be an integer. However, m is 
just a number parameterizing the form of the scalar field potential when we start with the action (|3.1I) obtained after 
dimensional reduction. We therefore find that m can be any real positive number in this context. The positivity 
needs to be assumed to keep the appropriate signature of the kinetic term for the scalar field (f>, or equivalently to 
keep the relation between <p and ip to be real. (Strictly speaking, the case with m < — n = —3 is also allowed.) Then, 
q can be regarded as a continuous parameter with its range < q < V3. 
The cosmological time r is related to 77 via 



dr = eS m+ ^ a / 3 dr]. 



(4.7) 



Recall that q^, the metric induced on the brane in the five[= n + 2]-dimcnsional model (|3.1(l . is related to by 
Eq. (|3.10l) . Hence the scale factor associated with the metric q^ is given by a — e mQ / 3 +/ 3 j an d therefore the Hubble 
parameter on the brane H := a~ 1 da/dr is given by H = (ma73 + /3')e" mQ/3 ~ Q . Substituting the above solution lETo^) 
into these expressions, we obtain 

2 ? 2 q(q 2 -3) 

a = {sin[(m- l)rj\} "C^-D {tan[(m- 1)77/2]} "(a*- 1 ) , 



H 



q(q 2 - 3 + 2gcos[(m - 1)77]) 
3(<Z 2 -3) 



{sin[(m — 1)77]} 



9(^-D {tan[(m - 1)77/2]} 
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FIG. 2: The motion of the scalar field in the potential. The behavior of the scale factor in the Jordan frame is also presented. 
The potential is positive for < q < 1 (left figure), whereas it is negative for 1 < q < \/3 (right figure). 

The relation between the coordinate time 77 and the cosmological time r (|4.7|l is not so obvious, but the asymptotic 

9(9 + 9) 

behavior can be easily studied. When 77 — > 0, we have r oc rj 9( -"+ 1 '> — > and 



a oc 779(9+1) oc t' , 

9(9+9) 9+9 
H oc r\ 9 (<j+ 1 ) oc a <j+ 3 . 



oc 77 



" 3(9+1) 



where the exponent p_ is defined below in Eq. (|4.8|l . For 77 := 77 — 7r/(m — 1) — > 0, the cosmological time is (locally) 



. 9(9- 



expressed as r oc 77 9 (9-i) . Therefore the range of the proper time r is infinite for the parameter region < q < 1 , 
while it is finite for 1 < q < \/3. In this limit 77 — 0, the scale factor, the Hubble parameter, and the scalar field 



behave like 



9(9-3) f t p+, (0<q<l), 
a oc 77 9(9-1) oc < , „ . _ 

\(T end -T)P+, {Kq<V3), 



9(9-9) 9-9 
H oc 77 9 (i-i) oc a 1- 3 , 



e oc 77 3 <i-«) 
In the above expressions, we have used 



+00, (0 < q < 1), 
0, (l<q<V3). 



P± : = , (4.8) 

4m + 9 ± ^3777(777, + 2) 

The ranges of p + and p_ are 1/(4 + y/S) < p+ < 1/3 and 1/3 < p_ < 1/(4 - y/S). 

The behavior of this solution is easily understood from the viewpoint of the four-dimensional effective theory 
described by the action (|3.13l) . as was discussed in Ref. (ijj. The potential (|3.14() with A& = and K = 1 is shown 
in FIG. [3 For 0<(7<l(0<m<l) the potential is positive, while for q > 1 (to > 1) the potential is negative. 
In the former case, the scalar field a starts at a — 00, climbs up the slope of the potential, turns around somewhere, 
and finally goes back to a — 00. In the latter case, a starts to roll down from a = 00. The universe expands for a 
period of time and eventually it starts to contract. Finally a falls into the bottomless pit within a finite time, where 
the universe ends up with a singularity. 

Here we note that the above picture based on the four-dimensional effective theory describes the dynamics in the 
conformally transformed frame in which the metric is given by q^, whereas we suppose that the "physical" metric is 
given by the induced metric on the brane q^. In principle, the cosmic expansion law can look very different depending 
on the frame we choose. The dynamics in the "physical" frame therefore can be very different apparently. However, 
the above discussion is still useful since the conformal rescaling does not change the causal structure of the spacetime. 
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B. Kasner type solutions with a cosmological constant 

The next example is a generalization of the Kasner- type spacetime including a cosmological constant Af,. Let us 
assume that the metric is in the form of 

g M dx A dx® = -dt 2 + e^S^dy^dy" + e 2m S ij dx i dx j , (4.9) 

where i and j again run from 1 to n, but here the metric of m-dimensional space is chosen to be flat (K = 0) because 
otherwise the solution with A& 7^ is not obtained analytically. The D[= n + 1 + m] -dimensional vacuum Einstein 
equations with a cosmological constant reduce to 

m (q + a 2 ) + n (S + /3 2 ) = (to + n)H$, (4.10) 
a + a (ma + n/3^j = (to + n)H$, (4-11) 
(3 + /3 (ma + n$) = (m + n)Hg. (4.12) 



From this we obtain two types of solutions (see Appendix B). One is a trivial solution, namely, /^-dimensional de 
Sitter spacetime, 

a = P = H a t. 

Of the other type is the following two solutions, 

e m a +nf3 = sinh [( m + n )H t], 



e 



, . 771 + 72 

tanh Hot 



±i/ ? 

(4.13) 



where q is the one that has been introduced in Eq. 14.4JI . The range of the time coordinate t is (—00,00) for the 
former de Sitter solution and [0, 00) for the latter non-trivial solutions. Applying the method discussed in Sec. Ill Al 
to these solutions, one can construct background solutions for a (D + l)-dimensional braneworld model. 

First we briefly mention the relation to the bulk inflaton model recently proposed by Koyama and Takahashi [25ll2^| . 
Identifying their model parameters A and S as 

9 , 2(to + 4) „ to + 4 l 2 H 2 

A = 46 2 - 8 /3 = — - — — , S=——— 5^r, 4.14 

1 to + 3 4(to + 3) 1 + £ 2 H 2 ' v ' 

we will find that our model is equivalent to theirs. The parameter to is supposed to take any positive number. Thus 
it follows that A varies in the same region, —8/3 < A < —2, considered in [25ll2^| . The background metric obtained 
by substituting the simplest solution a = (3 = H^t is indeed the case discussed in their paper. 




FIG. 3: The sketch of the behavior of the scale factor a (dashed lines) and the Hubble parameter H (solid lines) in the Jordan 
frame as functions of the proper time on the brane r. The left figure shows the solution of (14.611 with q > 1. The right figure 
describes the solution of (14.1311 with the plus sign in the exponent. 
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FIG. 4: The motion of the scalar field in the potential and the behavior of the scale factor in the Jordan frame. The solution 
found by Koyama and Takahashi pfll l2?SI | is described in the left figure. One of the non-trivial solutions with its behavior at 
the starting point a oc t p + is also shown in the left figure, while the other one with the exponent p_ behaves as shown in the 
right figure. 

Next we consider the cosmic expansion law. We start with the simplest case a = j3 = Hot. The dimensionally 
reduced metric (on the brane) is 

q &h dx & dx i = e 2mHot / 3 (-dt 2 + e 2H ° t 5 l:j dx t dx : >) . 

Introducing the cosmological time r and the conformal time r\ defined by 

dr = adr/ = e mH ° t/3 dt, 

the scale factor on the brane a = e ( m + 3 ) ff °*/ 3 j s written in terms of r or 77 as 

a cx T (m+3)/m cx v - {m+3)/3 . (4.15) 

Since 1 < (m + 3)/m < 00, power-law inflation with any exponent can be realized. 

Furthermore, we have non-trivial solutions i|4.13[l . The behavior of the solutions is as follows. At early times (t ~ 0), 

the scale factor behaves like a = e ma / 3 +' 3 ^ t 1 / 3 , and the cosmological time is given by dr ~ ^( m± -\/ 3m ( m + 2 )) / 3 ("i+ 3 ) ^ 
(and so r — > as t — > 0). Therefore, we have 

a~r p± , (4.16) 

withp-t introduced previously, which implies that the universe is not accelerated at early times. At late times (t — ► 00), 
we see that a — ► H t and (3 — > H t, and the solution shows power-law expansion as is given in Eq. I|4.15|l . 

A rather intuitive interpretation of the behavior of these three solutions can be made from the four-dimensional 
point of view again. The situation is summarized in FIG. 01 This time the potential is always positive. For one of the 
non-trivial solutions with the exponent p+ in (|4. 1(>|> . the scalar field starts to roll down from a — —00. For the other 
non-trivial solution with the exponent p_, the field starts to climb up the slope of the potential from a — +00, turns 
around somewhere, and rolls down back to a = +00. Suppose that the field a is increasing. Let us trace the evolution 
of a backward in time. If the kinetic energy is larger than a certain critical value, a will not have a turning point in 
the past. In this case a continues to decrease, reaching — 00. This corresponds to the case with the exponent p + . If 
the kinetic energy is lower, a will have a turning point. Then, we will have a — * +00 at t — > —00. This corresponds 
to the case with the exponent p_. The case of the power-law inflation l|4.15|l is, in fact, the marginal case between 
these two. In this case, a does not turn around. Therefore the evolution of a is similar to the case with p + . In any 
case, information about the initial velocity is lost as the universe expand. Therefore the late time behavior of the 
solutions is unique, and is given by Eq. I|4.15|l . 

V. COSMOLOGICAL PERTURBATIONS 

We consider cosmological perturbations in the (n + 2)-dimensional bulk inflaton models defined by l|2.12|l and l|2.13|l . 
The analysis of perturbations is complicated if we work in the original (n + 2)-dimensional models with a bulk scalar 
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field. Our (n + 2)-dimensional system, however, is equivalent to the (D + 1)[= n + 2 + rn]-dimensional one defined 
by (|2.2[1 and (|2.3|) . We will show that the perturbation analysis becomes very simple and transparent in the (D + 1)- 
dimensional picture, in which we just need to consider pure gravity without any matter fields. 
We begin with the following form of background metric: 

G AB dX A dX B = e 2 ^ (dz 2 - dt 2 + e^j^dy^dy" + e 2/3( %cfeW) , (5.1) 

where latin and roman indices in the lower case, respectively, rum m and n-dimensional subspaces, and the warp 
factor is given by e^W = £H /smh{H z). We assume that a(t) and f3{t) are chosen so that g^ B be a solution of the 
(n + 1 + m)-dimensional vacuum Einstein equations with a cosmological constant, 

m(& + a 2 ) + n(/3 + /3 2 ) = {m + n)H$, (5.2) 

a + a (ma + n0 \ + K(m - l)e~ 2a = (m + h)Hq , (5.3) 

(3 + /3 (ma + n/3) = (m + n)H„. (5.4) 

The background solutions with 7 M „ = 5^ (K — 0), Hq ^ and with K = ±1, Hq = were discussed in the 
preceding section. In the following discussions, we include more general cases with K = ±1 and Hq ^ 0. Although 
the background solution cannot be obtained in an explicit form for such non-flat compactifications with a cosmological 
constant, we will find that general properties of perturbations can be explored to a great extent. 
We write the perturbed metric as 

(G AB +5G AB )dX A dX B 

= e 2w {(l + 2N)dz 2 + 2Adtdz - (1 + 2<P)dt 2 + e 2a (l + 2S) 1 ^dy^dy v 

+e 20 [(1 + 2^)5 ij dx i dx j + 2E ij dx i dx i + 2Bidx l dt + 2C l dx i dz\ } . (5.5) 

These perturbations are assumed to be homogeneous and isotropic in the directions of the m-dimensional compactified 
space spanned by the coordinates From the assumption of isotropy, mixed components such as SG^tdy^dt and 
SG^idy^dx 1 are set to zero. Concerning the metric perturbations of the compactified space, therefore only the overall 
volume perturbation S is considered. After reduction to (n + 2)-dimensions, S is to be interpreted as the scalar field 
perturbation. Here, we also assume that the dependence on the n-dimensional coordinates x % is given by e lkiX . 

Metric perturbations are decomposed into scalar, vector, and tensor components based on the behavior under the 
transformation of the n-dimensional spatial coordinates x % in the following manner: 

Bi = ^-b s + bY, VBY -0, 

ik 

k l EY = k l Ej 3 = 0, S^Efj = 0. (5.6) 

The quantities with a superscript S, V and T represent scalar, vector and tensor perturbations, respectively. The 
perturbations SGab obey the linearized Einstein equations supplemented by boundary conditions at the position of 
the brane, 8K. A B \ z=Zb = where JCab is the extrinsic curvature of the brane. From the (D + l)-dimensional point of 
view matter sources are absent on the D-dimensional brane, and this makes boundary conditions considerably simple. 
Each component of the Einstein equations is presented in Appendix A. 

Here we would like to discuss the number of physical degrees of freedom in scalar, vector, and tensor perturbations. 
The transverse traceless tensor Ef- has (n + l)(n — 2)/2 independent components, each of which obeys a second 
order differential equation. For vector perturbations there are three variables BY , Ej and CY . The coordinate 
transformation has one vector mode, and correspondingly there is one vector constraint equation. Therefore we 
have only 1[= 3 — 1 x 2] vector remaining as a physical mode. Since a transverse vector has (n — 1) independent 
components, we find that there are 1 x (n — 1) degrees of freedom in vector perturbations, corresponding to the 
"graviphoton" . For scalar perturbations there are 8 variables, and the coordinate transformation has 3 independent 
modes. Since there are the same number of constraint equations, the number of physical modes is 2[= 8 — 3 x 2]. 
One of them corresponds to the bulk scalar field and the other corresponds to the "graviscalar" . In total, there are 
1 + (n + 2)(n - l)/2 = 2 + (n - 1) + (n + l)(n - 2)/2 physical degrees of freedom. The first "1" on the left hand 
side corresponds to the bulk scalar and the other (n + 2)(n — l)/2 degrees of freedom to those of (n + 2)-dimensional 
gravitational waves. 
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A. Tensor perturbations 



Since tensor perturbations are gauge-invariant from the beginning, they are in general easy to analyze. The 
equations for tensor perturbations are read from the {i, j}-component of the Einstein equations (|A5|1 as 



CE 



0. 



(5.7) 



where we have defined a differential operator 

C := d 2 + (ma + n^j d t + e~ 2f3 k 2 -d 2 z -(m + n){d z uj)d z . 
The perturbed junction condition implies that boundary conditions are Neumann on the brane, 



d,E. 



V \z=z b 



= 0. 



(5.8) 



Since the perturbation equations are manifestly separable, we write Efj = x(t)ip(z)Y? (x ) where is a transverse, 
traceless tensor harmonics. Then \(t) and ip( z ) obey 



X + (ma + nl3j x + (e 2f) k 2 + M 2 ) x = 0, 
d z ip + (m + n)(d z cj)d z i/j + M 2 ip = 0. 



(5.9) 
(5.10) 



Here M 2 is a separation constant and represents the squared Kaluza-Klein mass for observers on the D-dimcnsional 
brane. 

Now we discuss the mode function in the ^-direction ip( z )- Using a canonical variable ip :— e^ip, with 



m + n 3 



Eq. 15 . 1 Of) is rewritten into a Schrodinger-type equation, 

- d 2 ip + V(z)^j = M 2 ijj, 

where the potential is 



V(z) = ^ + 1)—^- 



sinh (Hqz) 



^H 2 -^S(z- Zb ). 



(5.11) 



(5.12) 



(5.13) 



The delta-function term is introduced so that ip automatically satisfies the boundary condition ^'(^fe) = 0. The 
presence of the zero mode, for which ip is constant in z, is obvious from Eq. I|5.10(l . From the asymptotic value of 
the potential V(oo) = /j. 2 Hq, we can say that there is a mass gap 5M — [iHo between the zero mode and the KK 
continuum. 

The z-dependence of the massive modes are given in terms of the associated Legendre functions by 



iM*) = c T (sinh(i? z)) 1 /2+M 



where c T is a normalization constant and 



(?-^(co S h(if oZt ))P;^(cosh(if o2 )) 



-Pyl^i^HHoz^Qj^licosHHoz)) 



' M 2 



(5.14) 



(5.15) 



These general properties of the mass spectrum and the mode functions in the z-direction hold irrespective of the 
specific form of the background solution a and /?. 

Let us move on to the time dependence of tensor perturbations. Using the cosmological time on the brane defined 
by r = J e ma l n dt, Eq. (|5.9|) is rewritten as 
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where one must recall that a = e ma / n +0 and H = a^ 1 da/dr. For the zero mode (M 2 — 0), this reduces to the 
equation for the tensor perturbations in the scalar-tensor theory defined by the action (|3.8|l . An apparent difference 
from Einstein gravity is the presence of the term (m/n)(da/dT). The Kaluza-Klein mass with respect to observers 
on the brane is expressed as 



m 2 KK {t)^e- 2ma{t ^ n M 2 , 

and so mkk = e~ ma l n ^jlHq for the lightest one, whereas the Hubble parameter at that time is given by 

H = e- ma / n (^-a + p) . 



(5.17) 



For a = = H t, this implies H = 2n 1 e 



-ma/n 



the other hand, when the background is given by Eq. I|4.13|) . we have H = 2n 1 e 



/j,H and therefore the mass gap and H are of the same order. On 



-ma/n 



/ii?o coth[(m + n)H t] and the 



mass gap can be very small compared to H , but only for a short period near t = 0. 

Despite its rather simple form, Eq. 15.9fl cannot be solved analytically in general. One exception is the case 
a = (i = H^t discussed in j^. In this case, Eq. I|5.9|l reads 

X + (m + n)H Q x + (e- 2Hot k 2 + M 2 ) X 

which, using the conformal time r\ = —e~ Hot /Ho, can be rewritten as 

drj 2 



1 



■ — n d 



V 



M 2 



This indeed has analytic solutions, 



drj 

Xocxi-nrHpi-kr)), 



= 0, 



x = o. 



(5.18) 

(5.19) 
(5.20) 



This is not a surprise because the background of the current model is just an AdS m +2+?i bulk with a de Sitter brane. 



B. Vector perturbations 

Next we consider vector perturbations. From perturbed junction conditions 5K. i |*=« 6 = and 8K, i t \ z=Zb — 0, we 
have 

d z EY\ =0, Cfl =0, d 7 B v A =0. (5.21) 

^ % \ Z — Zu 1 1 Z = Zu 7 ^ % \ Z — Zu v ' 



EY-Er + k^, BY^BY-^, CY = CY~d z g. (5.22) 



Under a vector gauge transformation x l — » x l = x L + , metric variables transform as 

r 

Thus we are allowed to set EY = by choosing an appropriate gauge. We expand the remaining variables by using 
the transverse vector harmonics as 

BY = BYY, CY = CY^. (5.23) 

For convenience, we introduce 

n := r 2 e m£,+ (" + ^e' m+ "> u (d z B - c) . (5.24) 
Then Eqs. (fASft . and IjAllj) are written as 

B = e - ma - n /} e -(.Tn+n,)u, Bg Q i ( 5 25 ) 

c _ e - ma +nf3 e -(m+n)w^ (5.26) 

It is easy to see that the remaining third equation is automatically satisfied if the above two equations hold. Substi- 
tuting these two into Eq. (|5.24JI . we obtain a master equation 

Q - (ma + n/3j Q + er 2p k 2 Q - [d 2 Q - (m + ri){d z u)d z Q\ = 0. (5.27) 
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This equation looks similar to the equation for tensor perturbations (|5.7(l . The difference is that the signatures of the 
terms containing first-derivatives such as f2 and £1' are reversed. From Eqs. ^5.2ip . the boundary condition for f2 on 
the branc turns out to be Dirichlet, 



ft|z=z 6 = 0. 



(5.28) 



Since the master equation l|5.27[l is separable, we write d(t, z) — x(t)£!(z). The canonical variable, Cl(z) :— e ^fl, 
again obeys a Schrodinger-type equation, 



d 2 Q + V(z)Q = M 2 Cl, 



with the potential 



Hi 



(5.29) 



(5.30) 



The crucial difference from tensor perturbations is the absence of the delta-function potential well. Because of this, 
there is no zero mode and only the massive modes with M 2 > V^oo) = ^ 2 Hq exist. The ^-dependence of the mode 
functions are given by 



n M (z) = c v (sinh{Hoz))V 2 -n 



Q-V2^( cosh ( ff o^))P_ 1 1 / / 2 ^( cosh (^o^) 



-P^Z( cosh ( H o z b))Q^ll(cosh(H a z^ 



(5.31) 



where c v is a normalization constant. When a = /3 = Hot, we can find an analytic solution for the time-dependence 
of the mode functions, which, using the conformal time, is given by 



C. Scalar perturbations 

1. Gauge choice, the boundary condition, and the mode decomposition 



Since scalar perturbations are more complicated, we begin with fixing the gauge appropriately in order to simplify 
the perturbed Einstein equations. We impose the Gaussian- normal gauge conditions 



N = A = C = 0. 



(5.32) 



Different from the case of vector perturbations, these conditions do not fix the gauge completely. In the case of 
scalar perturbations we need to take care of perturbations of the brane location. Here we make use of the remaining 
gauge degrees of freedom to keep the brane location unperturbed at z = Zb- In the Gaussian-normal gauge, boundary 
conditions on the brane for all remaining variables become Neumann: 



d z E\ 



d z S\ z 



d z B\ 



0. 



(5.33) 



Three of eight scalar perturbation equations are the constraint equations, and the other five are the evolution 
equations. First, let us examine the constraint equations IJA7I) . i|A13() . and l|A14l) . Eq. (| A 13(1 reduces to d 2 (& + n\l/ + 
mS) + d z Lud z {<$> + + mS) = 0. Taking into account the boundary conditions, this equation is once integrated to 
give 



d z ($ + nV + mS) = 0. 



Eqs. (|A7|) and l|A14|l reduce to 

d z 

o z 



kB + 2$ + 2 (ma + n/3) $ - 2maS - 2n/3* 
B + (ma + n$ + 2/3) B + 2e~ 2/3 fc 



= 0, 



--1)E 

n 



= 0. 



(5.34) 

(5.35) 
(5.36) 
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With the aid of Eq. (|5.34() . we find that all the perturbation equations are separable. Furthermore, the ^-dependent 
parts of these equations are the same as those of tensor perturbations with the same type of boundary conditions. 
Therefore we can expand all variables by using the same mode functions in the z-direction as those for tensor 
perturbations: 

E = E Q {t)^ {z) + Y J E M {t)i>M{z), $ = * (i)*o(z) + ^i(#M(4 (5.37) 

where i/^o is constant, ipu is given by Eq. I|5.14l) . and M 2 > ^ 2 Hq. Consequently, Eqs. (|5.34() . (|5.35() . and (|5.36() are 
automatically satisfied for the zero mode. For the massive modes these constraint equations give 

$ + + mS = 0, (5.38) 
kB + 2$ + 2 (ma + n$) $ - 2maS - 2n$^ = 0, (5.39) 



B + ( ma + n(3 + 2/3 ) B + 2e~ 2f3 k 



<]' + I - - i : / 

n 



= 0, (5.40) 



where the subscript M was abbreviated. Note that these three equations are nothing but components of the divergence 
of the metric perturbations, 

V A 5G zA = V A 5G tA = V A 8G lA = 0. 

In other words, the transverse traceless conditions are automatically satisfied if one impose the Gaussian normal gauge 
conditions except for the contribution coming from the zero mode. (V a 5G z a gives the traceless condition.) Below 
we discuss the KK modes and the zero mode separately. 



2. KK modes 

By using the constraint equations (|5.38() - (|5.40() . the Einstein equations (|A2|) . (|A3|) . IjAlOf) . (|A12|) . and (|A15|) are 

simplified to give 



= - 



- 1 kf3B + 2/3$ + 2(m + n)H^, 



CE = 2kf3B 1 

CS = kaB + 2a® + 2er 2a K(m - 1)(S - $) + 2(m + n)i/ 2 $, 
= -2 



[met + n(3 \ $ + 2k(3B + 2maS + 2nf3^ + 2(m + n)H^, 



CB = - I m'a + n(3 \ B — A(3B - Z B - 2(m + n)H'^B - Ae'^kfi® 



(5.41) 

(5.42) 
(5.43) 
(5.44) 

(5.45) 



where C is defined in Eq. I|5.8|l . Two of them give independent master equations for the massive modes, and the 
remaining three equations do not give any new conditions. With the aid of the constraint equations H5.38fl - H5.40fl . 
Eqs. (|5.41l) and l|5.44|l can be rewritten as 



-2(n- 2)P{a- (3) * + -/?$ 



-(3 



(m + l)d + n/3 — (m + n )H 2 \§ 



= -4/3$-2|(m+l)d + n/3 + 2/3 (m + l)d + n/3 -(™ + n)# 2 j$ 



-2n 



a - + 2$ (a - /?) 



(5.46) 



(5.47) 



Unfortunately, except for the simplest case (to be discussed later) we do not know how to disentangle these two 
equations, although there is no problem in solving these equations numerically. Once we could solve these coupled 
equations, the other variables S, B, E are easily determined just by using the constraint equations. 
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3. Zero mode 



To discuss the zero mode, it is useful to look at the cosmological perturbations in the corresponding (n + 1)- 
dimensional theory defined by the action i|3.13[l . in the case of the (n + l)-dimensional Friedmann universe with 
a single scalar field, there is only one physical degree of freedom in scalar perturbations. One can derive a second 
order differential equation for one master variable [44J- Back in the braneworld context, the background metric and 
its zero-mode perturbations are also described by the same effective action (|3.13|) . Therefore, the analysis of the 
zero mode is no different from the conventional [n + l)-dimensional cosmological perturbation theory. Below we will 
explain this fact more explicitly. 

To begin with, we consider (n + 1 + m)-dimensional spacetime whose metric is given by 



ds 2 = -(1 + 2<$>)dt 2 + e 2a (l + 2S)~j llv dy >l dy v + e 20 [(1 + 2V)5 ij dx i dx j + 2E ij dx i dx j 



2B t dx l dt] 



(5.48) 



where only scalar perturbations are imposed and they are again assumed to be homogeneous and isotropic with respect 
to the m-dimensional compactified space spanned by y 11 . Then, the perturbed Einstein equations ^R^g = (tn + 
n)H 2 5g AB become identical to Eqs. <|A2I) . (|A3|) . (jA10|) . (|A12|) . and (|A15|I with N, A, C, and the terms differentiated 
by z dropped. Hence, it is manifest that the analysis of zero-mode perturbations in our (n + 2 + ni)-dimensional 
spacetime is equivalent to that of the above system. 

As for perturbations in (n + 2 + rn)-dimensional spacetime, we have already fixed the gauge by imposing three gauge 
conditions H5.32[l . However, these gauge conditions do not fix the gauge completely. As is manifest from Eqs. I|A17|I . 
gauge transformations satisfying £ z = and £ s = £* = do not disturb the conditions Ij5.32|) . On the other hand, 
on the (n + l)-dimensional side there are two scalar gauge transformations 

x i -s- x 4 = x* + k l £ s /ik. 

The transformation of metric variables under these gauge transformations is the same as that obtained by setting 
£ z = and £ s = = in the last five equations in (|A17J| . 

If we think of the size of compactified dimension S as a scalar field in (n + l)-dimensional spacetime, the system 
reduces to a conventional (n+ l)-dimensional model with a scalar field. In the conventional cosmological perturbation 
theory, $ and 'J in the longitudinal gauge (B = E = 0) are known to be convenient variables. Here one remark is 
that we need to take account of a conformal transformation to map the theory to the conventional (n + l)-dimensional 
one, 



dS 2 



^ (a+S ) . e ^( a +S) ds 2 = e ^ {a+S ) ds 2 



(5.49) 



which follows from the discussion in Sec. IIIII Then, the variables corresponding to the so-called Sasaki-Mukhanov 
variables are 



$ := $ 
* := * 



n — 1 

m 

n — 1 



S, 

s. 



in the longitudinal gauge. 

Eliminating N, A, B, C, E and the terms differentiated by z, Eq. 



2)* = 0. 



Similarly, from Eqs. l|A10j) . (|A2|> and l|A15|l . we have 

(n - 1)* + (n - 2) ma + (n - 1)0 



3l becomes 



1 



n- 1 



-maS. 



where £ f ■= I 



= -2(n - l)/3* + 2(m + n)H$$, 
C a S = -2(n - l)a* + 2(m + n)Hfo - 2K(m - l) e - 2a {<5> - S), 

n/j) / + e~ 2f3 k 2 f. Combining all these, we obtain the EOM for 4": 



(5.50) 
(5.51) 

(5.52) 



(5.53) 



A,* - 2 ( /3 + | J 4 + 2(7i -2) ( (3-^)^ = 0. 



(5.54) 
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Using the conformal time r\ = J e @dt, we can rewrite this into a more familiar form as 



(n - 1)H - 2— 
a 



+ fc 2 * + 2(n - 2) ( H' - H— ) * = 0, 



(5.55) 



where we have denned 



H := 



(In a)' = [ma + (n 

n — 1 



(5.56) 



and 



*/(™ i s the scale factor in the Einstein frame. 



Since there is a mass gap between the zero mode and the massive modes in general in our models except for a short 
period in the cases of K ^ 0, the massive modes would not be excited easily. Hence, the behavior of the zero mode is 
especially important. Since we found that the zero mode can be described by the corresponding (n + l)-dimensional 
conventional cosmology, it can be easily analyzed in general. 



4- Exactly solvable case 

Let us consider the simplest background given by a = (3 = H$t with K = 0. In this special case, scalar perturbations 
including the KK modes are solved exactly. The most remarkable advantage of our approach may be that z-dependence 
of the modes can be derived for a general background as we did in the earlier part of this section. The time dependent 
part, which is usually non-trivial especially for the KK modes, is also solved easily as shown below when a = ft = H^t. 

Substituting a — f3 — H t into Eq. (|5.47|l . the equation for $ is decoupled first, 

£$ = -4iJ $ - 2(m + n + 2)H$$. (5.57) 

By assuming ^-dependence given in (|5.14|) . we expand as $ = ^m^m- Then, using the conformal time, the above 
equation is rewritten as 



d 2 m + n + 3 d ,„ 1 f M 2 . 

+ k + T T7T + 2 ( m + n + 2 ) 



-r-7 r " - i 77T -w« i , $m = 0. (5.58) 

a?7 r\ ax] r] z \ Hq j 

The solution is given in terms of the Hunkel function by 

$M = cf(-r?)V^ (5-59) 

with 

P = {-riYH\l\-kri), 

where cf is a constant and fi and v were denned in Eqs. (|5.11|l and l|5.15|l . Then, substituting this into Eq. (|5.4U|) 
with the aid of Eq. (|5.39(l , Bm is immediately obtained as 

B M = -24k- 1 H [{-r,fp' + (2fx - p] ■ (5.60) 

The result is consistent with the evolution equation for B (I5.45|l . 

Eqs. I|5.42|l . (|5.43|) and l|5.45|l are combined to give a simple equation, 

C[^m + E M In - S M ] = 0. 

The operator C is the one that appeared in tensor perturbations, and so the mode solutions were already known: 

* M + E M jn - S M - cf p, (5.61) 

where cf is another constant. Substituting 4>m and Bm, the constraints l|5.38l) and 15.40fl reduce to two algebraic 
equations for ^m, Em, and Sm as 

n^M + mSM = -cfi-n) 2 p, (5.62) 
*m + (1/n - l)E M - -cf {-nfp - cf fc- 2 [(2/i - l)r,p' + {v 2 - 3p 2 + 2p)p] . (5.63) 
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Solving Eqs. (|5.61(l . I|5.62[l and (|5.63|l . we obtain the expressions for 'I'm, Em, and Sm- Thus, all the metric variables 
can be analytically solved. Note that one of the above two independent solutions were already obtained in Ref. [25ll26l| . 
The zero-mode solution is also easily obtained. In this background, the master equation becomes 

m + n ~H ' + fc 2 ^ = 0, (5.64) 
V 

and the solution is 

t< = c s (- V r- 1 Hl 1 l 1 (-kr ] ). (5.65) 



VI. CONCLUSIONS 



We have shown that a wide class of braneworld models with bulk scalar fields can be constructed by dimensional 
reduction from a higher dimensional extension of the Randall-Sundrum model with an empty bulk. The sizes of 
compactified dimensions translate into scalar fields with exponential potentials both in the bulk and on the brane. 
We have mainly concentrated on models with a single scalar field, which include the power-law inflation solution of 

Ref. mm. 

First we have investigated the evolution of five[= n + 2]-dimensional background cosmologies, giving an intuitive 
interpretation based on the four[= n + l]-dimensional effective description. Then we have studied cosmological 
perturbations in such braneworld models. Lifting the models to (5 + m)[= n + 2 + m]-dimensions is a powerful 
technique for this purpose. The degrees of freedom of a bulk scalar field in (n + 2)-dimensions are deduced from a 
purely gravitational theory in the (n + 2 + m)-dimensional Randall-Sundrum braneworld, which consists of a vacuum 
brane and an empty bulk. We would like to emphasize that the analysis is greatly simplified thanks to the absence 
of matter fields. From the (n + 2 + m)-dimensional perspective, we have derived master equations for all types of 
perturbations. We have shown that the mode decomposition is possible for all models which are constructed by using 
this dimensional reduction technique. Moreover, the dependence in the direction of the extra dimension perpendicular 
to the brane can always be solved analytically. 

As for scalar perturbations, there are two physical degrees of freedom for the massive modes and the equations are 
not decoupled in general. For the zero mode, however, the situation is equivalent to the standard four-dimensional 
inflation driven by a single scalar field. Hence, only one degree of freedom is physical. Therefore we end up with a 
single master equation. To sum up, our "embedding and reduction" approach enables systematic study of cosmological 
perturbations in a class of braneworld models with bulk scalar fields. 

In this paper, we have not discussed quantum mechanical aspects. In order to evaluate the amplitude of the 
quantum fluctuations, the overall normalization factor of the perturbations must be determined. For this purpose, 
one needs to write down the perturbed action up to the second order written solely in terms of physical degrees of 
freedom as is done in the standard cosmological perturbation theory. We would like to return to this problem in a 
future publication. 

In this paper our investigation is restricted to the parameter region m > 0, where m is the number of compactified 
dimensions. If m > 0, a singularity could exist at z = oo, but it is null. For m < — 3, we have the right sign for 
the kinetic term of the scalar field and so it is possible to consider such models. In this parameter region, however, 
there is a timelike singularity at z — and therefore we need a regulator brane to hide it. This case includes the 
cosmological solution of heterotic M theo ry Il3 | (which corresponds to m = —18/5), and the analysis of cosmological 
perturbations in such two-branes models 17] would be also meaningful. This issue is also left for future work. 
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APPENDIX A: PERTURBED EINSTEIN EQUATIONS AND GAUGE TRANSFORMATIONS 



In this appendix, we write down the components of the perturbed Einstein equations 

C 1 ^ t- i„ „ 1 + m + n 



J) 



±n6G AB - \v A V B 8G 



P 



SG 



AB- 



The perturbed quantities are decomposed into scalar, vector, and tensor components whose basic definitions are given 
by Eq. (|5.6[1 . Note that in the following expressions no gauge conditions have not been imposed yet. 



{i, j}-component 

1 r 



hij + (ma + n0)hij - - (to + - e~ 20 d k d k h li + 2e~ 2fs d k d {l h j)k 



d(iBj) + (ma + n/3)d(iBj) — d^C'.^ — (m + n)u>'d^Cj) 
+%|2(m + n)H 2 (N - $) + 2(1 + m + n)uo" N + lj'N' 



8a (p d k B k - uo'd k C k 



ma + (to + 2n)/3 



J A - (3 A' - J A 



/ ($ + mS + n*)' + /3 (N - $ + mS + - e" 2 ^^ (AT + $ + toS + = 0, 



(Al) 



where hij = 2\&&y + 2Eij and the dot (prime) denotes d/dt (d/dz). (We use prime to denote differentiation 
with respect to z only in Appendix A.) 

Trace Part: 



* + (ma + n/3)^ + e 2f3 k 2 ^ - - (m + n)uj'^>' 



*+ --1 I f? 
n i \n 

+ ^2(m + n)H 2 (N - $) 
+2(1 + m + n)uj"N + uj'N' - ma + (m + 2n)(3 



B+ (ma + 2n(3) B - C' - (m + 2n)JC 



J A - /3 A' - J A 



-J (<$> + mS + n*)' + (3 (N - $ + mS + n^y) + e~ 2f3 k 2 (N + $ + mS + n*) - = 0. 

J n 



(A2) 



Trace- free Part: 



E + (ma + n(3)E + e - 2f3 k 2 E - E" - (m + n)u>'E' + 2e~ 2f3 k 2 









(=-')-] 







+k 



Vector: 



Tensor: 



{z, i}-component 



B + (ma + n/3)B~\ - k[C' + (m + n)u)'C] - e~ 2(3 k 2 (N + $ + mS + n#) = 0. 

EY + (ma + n0)Ej - Ej" - (m + n)JEj' 

+k \bY + (ma + n/3)BY] - k [Cf' + (m + n)w'CY ] = 0. 

d 2 + (ma + n[3)d t + e^k 2 - <9 2 - (to + n)Jd} Ej- = 0. 



(A3) 



(A4) 



(A5) 



(d j hij) - (diAy- [ma + n{3 - 2(3) d t A + d j d l C j - d j d j C l 
~2d t ($' + mS' + n^') + 2(m + n)uj'd l N 



+e 
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Ci — B 



?/ + (ma + n/3 + 2/3) (di 



B' 



0. 



(A6) 



Scalar: 



k 



2*' + 2 ^ - 1^ E' - A - (ma + n/3 - 2/?) A - 2 ($' + mS' + n*') 



Vector: 



{t, i}-component 



+2{m + n)uj'N 
kEY' + k 2 C'Y + e^ 



C-B' + (ma + n/3 + 2/3) ((5* - B') 
C'Y - BY' + (ma + n/3 + 2/j) (cf - B Vl 



0. 



= 0. 



(d i h ij y+ (d t A)' + (to + n^'^A + 0*5^- - ffdjBi 

+2 (md + n/3 - /?) <9 4 $ - 2d t (n + mS + ntf) + 2d t 0N + m (/? - d) S 



+e 
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-Bi" + C- + (to + n)J (di - Bi ,y j 



0. 



Scalar: 



2* + 2 Q - 1^ E + A' + (m + n)JA + 2 (ma + n/3 - /3) $ - 2 (/V + mS + n*) 



+2/3 + 2to (/3 - d) 5 



Vector: 



kE v + k 2 B t + e 2 ? 



B" + C' + (m + n)J (C - 5') 
BY" + CY' + {m + n)J (CY-BY 1 ) 



0. 



0. 



{t, i}-component 



k (b + 2(3B - u/c) - (ntf + mS)"+ w' ($ + n* + mS)' + A' + (1 + to + n)u/A 

+ (m + n)u/$' - e~ 2,3 fc 2 $ + (ma + n/3j $ - 2n/3* - 2mdS' 
-/V - u/7V' - 2(1 + to + n)u"N - 2(m + n)H$(N - $) = 0. 



{z, z}-component 



TV + e- 2/3 fc 2 iV + fc (C + - ($ + n* + mS)" - J ($ + n* + toS*)' 

- (i' + uj'Aj - (ma + n/3) + w'A - TV) + (1 + m + n) (w'JV' + 2J'N) = 0. 



{£, z}-component 



{fi, ^}-component 



k (C + 2(3C + B') + e- 2f3 k 2 A + 2 (ma + n/3) $' + 2u/(m + n)N 
-2 (n& + mS") - 2n/3*' - 2mdS" - 2(m + n)i7 2 A = 0. 



5 + (2md + n/3) S - S" - (2m + n)u/S' + [er 2fJ k 2 - 2 e - 2a K(m - 1)] S 

+nai> - d$ - u/ ($' + n*') + 2e 2a K(m - 1)$ + k (JC - aB) 
+aN + lu'N' + 2(1 + to + n)uj"N - J A - aA' - [(2m + n)a + n/3] J A 
+2(m + n)Hl(N - $) = 0. 
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Lastly we summarize the gauge transformations of the metric variables. Under a scalar gauge transformation, 

t -► + 

z -> z = z + £ 2 , 

a:* -» + ^rf, (A16) 

the metric variables transform as 



A - 


N = 


N- 




A - 


■ A = 






C - 


■ C = 




-e-*Pk?-Z 8 ' 


B - 


■ B = 


B - 




$ - 


■ $ = 


$ - 




* - 


• § = 


* - 


n 


E - 


E = 


E-\ 




S - 


5* = 


S- 





APPENDIX B: MULTI SCALAR FIELD GENERALIZATION 



Let us give the generalization of the Kasner type metric discussed in Sec. IIVI First we generalize the case without 
Ab but including the curvature for one of the spatial sections: 

M 

g MN dx M dx N = e 2 ^) [_ d?? 2 + 7/B , dlW ] + e 2 ^)5%dx M dx N , 

i=l 

where 5^ N is the metric of a j.;-dimensional flat space and is the metric of a m-dimensional maximally symmetric 
space. As before Here we identified j% with n. If we compactify {fh + ^^L^Ji) dimensions leaving (n+ 1) dimensions, 
the compactified space is divided into N sectors having different scale factors. Then the resulting cosmology after 
dimensional reduction will possess N scalar fields. 
The set of vacuum Einstein equations becomes 

e 2a RJ = Jilf ~ a'u + u' + ma" = 0, (Bl) 

e 2a R^ = <y \K(m - 1) + a'u + (to - l)a' 2 + a" = 0, (B2) 

e 2a R£S$ L = S$ L Hu + (to - 1)«' 7 .: + 7-'] - 0, (B3) 



where we have introduced 



From Eq. I|B3|) . we obtain 



: ( B4 ) 



u 2 + (to - l)a'u + v! = 0. (B5) 



Then, it is easy to see that Eqs. <|B2(1 and (|B4(1 are equivalent to Eqs. (|4.2|l and (|4.3|) in the example of Sec. IIV Al by 
identifying u with n(3' . Therefore the solution of Eqs. I|B2|I and l|B4|l for K = 1 is written as 

±(to — l)o , r , „ , o ,„ , 

a = cot[(TO - l)n] T . u - 7VT ' ( B6 ) 



sin [(to — 1)??] ' sin[(m — I)?/] 



Since these two equations <|B2ll and (|B4f) do not have dependence on the number of dimensions, q has not been fixed 
yet. Substituting this solution into Eq. l|B3j) . we obtain 

(m-l)^cot[(m-l)j7]+ 7 f = 0, 
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which is integrated to give 



7i 



(to - l)cj 
sin[(m — 1)77] 



where Cj is an integration constant. Substituting this into the remaining equation (|B1|I and the definition of u (|B4|l . 
we find that the solution is given by 

±9 



cot 



m — 1 



-v 



with 



to — 1 
tan [ 77 



V i,c? = — , V j,c 4 = ±g. 

to — 1 



The next is a generalization of the Kasner-type spacetime including a cosmological constant Let us assume that 
the metric is in the form of 



9mn 



dx M dx N = ~dt 2 + J2 e 2 ^8 { ^ N dx a dx b 



(B7) 



Here all the spatial sections are taken to be flat (Ki = 0), because otherwise an analytic solution with ^ cannot 
be found. The Einstein equations Rmn[9] — NH^gMN, with N := D — 1 = J2 n i reduce to 



7« + = nhI, 



(B8) 
(B9) 



where dot denotes differentiation with respect to t. These equations admit a trivial solution of £>-dimensional de 
Sitter spacetime, 

7l = Hot. (BIO) 
There is another type of non-trivial solutions. From Eq. (|B9|I we find that u := obeys 

ii + u 2 = N 2 Hq. 

The solution for this equation is 

u = NH coth(NH t). 

Substituting this into Eq. (|B9|) . we obtain 



[sinh(JVif t)7i] ' = NH 2 svah(NH t). 
This can be easily integrated and integration constants are determined from Eq. I|B8|) . Then we have 

. H [cosh(NH t) + a] 
7i ~~ sinh(iViJoi) ' 
5>c,=0, ^j iC 2 = 7V(7V-l). 



Finally, integrating Eq. (|B11() . we obtain 

e* 7 ' = smh.(NH t) 



tanh 



/NHo t 



(Bll) 
(B12) 

(B13) 



Integration constants were removed by rescaling the spatial coordinates. There are only these two types of solutions 
(|B"TU|) and (|B"T3|) for Eqs. JEHJ and JH5J. 
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